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Abstract 

In this paper, we investigate the two component long wave short wave resonance interaction 
(2CLSRI) equation and show that it admits the Painleve property. We then suitably exploit 
the recently developed truncated Painleve approach to generate exponentially localized solutions 
for the short wave components and 5(2) while the long wave L admits line soliton only. The 
exponentially localized solutions driving the short waves S^^^ and S^'^^ in the y direction are endowed 
with different energies (intensities) and are called "multimode dromions". We also observe that 
the multimode dromions suffer intramodal inelastic collision while the existence of a firewall across 
the modes prevents the switching of energy between the modes. 
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I. INTRODUCTION 



Recent investigations of the integrable coupled nonlinear Schrodinger equation, namely 
the celebrated Manakov model, and the observation of intensity redistribution in the collision 
of sohtons [1-5] have clearly pointed out their potential usage in the field of optical commu- 
nications and have virtually set in motion the process of designing an all optical computing 
machine. In particular, the vector solitons undergoing energy sharing collision identified in 
the coupled NLS equation turned out to be the key in the growing list of alternatives to 
the paradigm of soliton based chips, at least for specialized applications including quantum 
computing [6], DNA computing [7] and dynamics based computing based on chaos [8]. 

It is known that the identification of dromions [9,10] in the Davey-Stewartson equation 
which is a (2+1) dimensional generalization of the NLS equation has given the much needed 
impetus to the investigation of (2+1) dimensional integrable models. These dromions which 
are localized exponentially in all directions are essentially driven by certain lower dimensional 
arbitrary functions of space and time. In fact, such lower dimensional arbitrary functions of 
space and time have consolidated the concept of integrability of the associated dynamical 
systems in (2+1) dimensions besides being tailor made for the construction of various kinds 
of locahzed solutions. Reflecting on the flurry of activities taking place in the held of optical 
communication ever since the identification of shape changing collision of vector solitons 
in the coupled NLS equation and the rapid strides made in the field of (2+1) dimensional 
nonlinear partial differential equations (pdes) after the observation of dromions in the Davey- 
Stewartson I (DSI) equations, one would be tempted to look for the possibility of identifying 
the counterparts of vector solitons in (2+1) dimensions as well. 

In fact, the recent derivation of the two component long wave short wave resonance 
interaction (2CLSRI) equation in the context of the interaction of nonlinear dispersive waves 
on three channels [11] has only fuelled the anticipation to look for such localized excitations. 
This is also further supported by the study of collision behaviour of plane solitons admitted 
by the 2CLSRI equation recently [12]. In this paper, we investigate 2CLSRI equation and 
confirm its Painleve property. We then suitably employ the recently developed truncated 
Painleve approach [13-16] and generate multimode dromions. It should be mentioned that 
this is the first time the existence of exponentially localized solutions has been reported in 
a vector (2+1) dimensional nonlinear pde. Finally, we also study the unusual interaction of 



2 



multimode dromions. 

We now consider the two component long wave short wave resonance interaction 
(2CLSRI) equation in the following form 

i(5« + 4^))-5«+L5«=0, (la) 
^(5r + 5f)-5i+L5(^)=0, (lb) 

Lt^2i\S(X + \S^X)- (Ic) 

The above equation is the two component analogue of the long wave short wave resonance 
interaction equation investigated recently [13]. In eq.(l), S^^^ and S^'^^ represent short waves 
while L denotes a long wave. In particular, it explains the interaction of long interfacial wave 
(L) and a short surface wave (S) in a two layer fluid. This equation has been investigated 
recently and hue solitons have been generated [11,12]. 



II. SINGULARITY STRUCTURE ANALYSIS 

We now rewrite the above equation by putting S'*-^-' = p, S'*-^-'* = q, S^'^'^ = r, S^'^'^* = s as 

i{Pt + Py) - Pxx + Lp = 0, (2a) 

-i{Qt + Qy) - Qxx + Lq^O, (2b) 

iin + Ty) - r^x + i^r = 0, (2c) 

-i{st + Sy) - s^^ + Ls = 0, (2d) 

Lt^2{pq), + 2{rs)^. (2e) 

We now effect a local Laurent expansion of the variables p, g, r, s and L in the neighbourhood 
of a noncharacteristic singular manifold (j){x^ y, t) — 0, (p^ ^ 0, 4>y ^ 0. Assuming the leading 
order of the solutions of eq. (2) to have the following form 

p = por, q = ?o0^ r = ro(j>\ s = So0^ L = , (3) 

where Pq, Qq, Tq, Sq and Lq are analytic functions of {x, y, t) and a, /9, 7, 5 and m are integers 
to be determined, we now substitute (3) into (2) and balance the most dominant terms to 
obtain 

a = /3 = 7 = 5= -l,m = -2, (4) 
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with the condition 



Poqo + roso = Lo = 24)1 



(5) 



Now, considering the generahzed Laurent expansion of the solutions in the neighbourhood 
of the singular manifold 



q^qo4>^ + ... + qj(p^-^^ + ..., 
r = rgcfP + ... + TjcfP'^'^ + 

S = Sq4>^ + ... + Sj(tP+^ + 

L = Lo0'^ + ... + L,-<^'+'^ + ..., 



(6a) 
(6b) 
(6c) 
(6d) 
(6e) 



the resonances which are the powers at which arbitrary functions enter into (6) can be deter- 
mined by substituting (6) into (2). Vanishing of the coefficients of {(j>'~^,(f>'~^,(l>'~^,(j>>~^,(j>'~^) 
leads to the condition 
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From equation (7), one gets the resonance values as 

j = -l, 0,0,0, 2, 3, 3, 3, 4. 



(8) 



The resonance at j = -1 naturally represents the arbitrariness of the manifold (f){x, y, t) = 
0. In order to prove the existence of arbitrary functions at the other resonance values, we 
now substitute the full Laurent series 



3 3 3 

s = so4>^ + «J■</^'^^ L = Lor + E ^^"f^^^ 
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into equation (2). Now, collecting the coefficients of (cf) ^,0 ^) and solving 

the resultant equation, we obtain equation (5), implying the existence of a resonance at 
i = 0,0,0. 

Similarly, collecting the coefficients of and solving the resultant 



equations by using the Kruskal's ansatz, ^(x, y,t) = x + i^iy, t), we get 

Pi = ^[iPoi'tpt + ^y) - 2p0a;], (lOa) 

9i = ^ [-ko (A + i^y) -'^Qox], ( 10b) 

= ^[^Mi^t + i'y) - 2roa;], (lOc) 

Si = ^[-isoiA + i^y) - 2soa;], (lOd) 

Li = 0. (lOe) 

Collecting the coefficients of we have 

i{pot + Poy) - Poxx + L0P2 + Lipi + L2P0 = 0, (11a) 

-i{qot + qoy) - qoxx + ^0^2 + ^igi + -^^2^0 = 0, (lib) 

ii'^ot + roy) - roxx + Lor2 + LiVi + Laro = 0, (11c) 

-i{sot + Soy) - so^^ + LqS2 + I/iSi + L2S0 = 0, (lid) 

Lit = 2[poxgi + qixPo + Pixqo + Piqox] + '^[roxSi + rixSo + si^tq + nso^] = 0. (lie) 



From (11a), (lib), (11c) and (Hd), we can eliminate L2 to obtain the following three 
equations for the four unknowns p2, ^2, and S2, 

^o(Pog2 - qoP2) - {PoqOxx - qoPOxx) - i{Po{qot + qOy) + qoiPot + POy)) = 0, (iif) 

Lo{por2 - roP2) - {poroxx - ropo^x) - i{-Po{rot + roy) + ro(pot + Poy)) = 0, (llg) 
Lo{poS2 - S0P2) - {posoxx - soPoxx) " i{po{sot + Soy) + So{pot + POy)) = 0, (llh) 

which ensures that one of the functions p2, q2, T2 or S2 is arbitrary. Obviously L2 itself 
can be obtained from any one of the four equations (11a), (lib), (11c) or (lid). Similarly, 
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collecting the coefficients of we have 



KPU + P2A) + iiPly + P2V^j/) - {Plxx + 2p2a;) + -^2Pl + -^SPO = 0, (12a) 

-i{qit + g2^t) - + q2'>Py) " (gixi + 2g2x) + -t'2gi + -t/ago = 0, (12b) 

i{rit + r2il^t) + i(?^i2/ + ^2^y) - {rixx + 2r2a;) + -^2^1 + Lgro = 0, (12c) 

-i{su + S2V't) - K^ly + «2V'2/) - («lxa; + 2S2x) + L2S1 + L3S0 = 0, (I2d) 
L-2t + i^sV't = 2[poa;g2 + {pix + ^2)?! + (p2a; + ^3)^0 

+qoxP2 + (gix + q2)p\ + (g2x + 53)^0] + 

2[roa;S2 + (rix + r2)si + (r2a; + rz)sQ 

+Sox?'2 + (six + S2)ri + (s2x + S3)?'o]- (12e) 

Equations (12a), (12b), (12c) and (12d) can be solved for L3 as 

= —{-i{pn + P2i^t) - i{Ply + P2'4^y) + {pixx + 2p2a;) " -t'2Pl), (I2f) 
I'D 

^3 = — (^(9u + ?2'0«) + ^{Qiy + ?2V'2/) + {qixx + 2g2x) - -^^2^1), (12g) 

L3 = — (-«(r-it + r2^t) - i(rij^ + r2ijy) + {ri^x + '^r2x) - ^-2^1) (12h) 

L3 = — (i(sit + 52^'*) + ^(Sly + S2V'2/) + + 2*20;) " L2S1). (I2i) 

So 

Making use of eqns. (5), (10) and (11), we find that the right hand sides of eqs. (12f), 
(12g), (12h) and (12i) are equal. This implies that we are left with two equations for five 

unknowns. So, any three of the five coefficients pa, qs, r^, S3 or L3 are arbitrary. Now, 
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collecting the coefficients of {4>, 0, 0, 0), we have 

i{p2t + 2p3ll;t) + i{p2y + ^pstpy) - {p2xx + ^PZx + 6^4) 

+L0P4 + L2P2 + L-iPl + L4P0 = 0, (13a) 
-i(g2t + 2g3?/^t) - i(Q2j/ + '^qzi^y) - {(hxx + 4g3x + 6^4) 

+Log4 + -^^2^2 + -C'S?! + ^4^0 = 0, (13b) 

i(r2t + 2r3V^t) + i(r2j; + "Ir-ii^y) - {"^ixx + 4r3a; + 6r4) 

+Lor4 + L2r2 + L^^ri + L4ro = 0, (13c) 

-%{s2t + 253^^*) - i{s2y + 2s3V'j/) - {&2xx + 4^3^ + 6S4) 

+L0S4 + -^2^2 + -^aSi + ^4*0 = 0, (13d) 

L3t + 2L4^t = 2[poxg3 - P0?4 + {P\x + P2)g2 

+ (P2a; + 2p3)gi + (PSo; + 3^4)^0 + ^Oa^PS " 90^4 
+ (gia; + g2)p2 + {q^x + 2g3)Pl + (fe + 3g4)po] + 

2[roxS3 - '^oS4 + {rxx + ^2)^2 + (r2a; + 2r3)si + 
{r-Ax + 3r4)so + 50x^3 - s^r^ 

+(sia; + S2)r2 + {s2x + 2s3)ri + {s-^x + 3s4)ro] . (13e) 

By multiplying (13a) by (13b) by (13c) by Sq, (13d) by Tq and adding the resultant 
equation, we obtain an equation which is same as (13e). This means that we have only four 
determining equations for five unknowns. So, any one of the five functions p4, ^4, r4, S4 or 
L4 is arbitrary. One can proceed further to determine all other coefficients of the Laurent 
expansions (9) without the introduction of any movable critical manifold. Thus, the 2CLSRI 
equation indeed satisfies the Painleve property. 

III. TRUNCATED PAINLEVE APPROACH AND LOCALIZED SOLUTIONS 

To generate the solutions of 2CLSRI equation, we now suitably exploit the results of 
the leading order behaviour by employing the truncated Painleve approach. Truncating the 
Laurent series of the solutions of eq. (2) at the constant level term, one obtains following 
the Backlund transformation 
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Po . 
P^-r+Pi, 



+ ri, 



So 

4> 



(14) 



Assuming the following seed solution 



Pi = gi = ri = si = 0, L2 = L2{x, y), 



(15) 



we now substitute (14) with the above seed solution (15) into equations (2) and obtain 
(5) by collecting the coefficients of , 4>~^ , (f)~^ , (f)~^ , (f)~^) . Gathering the coefficients of 
(0~^, we have the following system of equations 



- ipo(l>t - ipo4>y + '2pox4>x + Po(f>xx + Lipo = 0, 

i%(l>t + «go0y + '^qOx4>x + Qoff^xx + LiQq = 0, 

-zro0t - iro(l)y + 2rQ^(j)^ + ro(/)j.a; + -^i^^o = 0, 
Lot - ^i*/"* = 2(po?o + ?'oSo)x- 



From equation (16e), we have 



Li = 2 



<t>t 



Using (17) in (16a-16d), the variables Pq, Qq, Tq and Sq can be solved as 

20a; <l>tx 



Po = Fi{y,t)exp 



ro = F2{y,t)exp 



So = F2{y,t)eiqp 



1 /• i(0t + (Py) + 0^ 

2 

1 f -K4>t + 02/) + 4>xx - 



^dx 



20a: 4't: 



-dx 



An 



+ 



^dX 



-dx 



(16a) 
(16b) 
(16c) 
(16d) 
(16e) 



(17) 



(18a) 



(18b) 



(18c) 



(18d) 



where Fi{y,t) and F2{y,t) are lower dimensional arbitrary functions of y and t. 
Substituting (18) in (5), we obtain the condition 

F,{t-yf^cl>t-F,{t-yY. (19) 

Again, collecting the coefficients of {4>~^, (f)~^, (f)~^, we have 

iPot + iPoy - Poxx + L2P0 = 0, (20a) 

-igot - iqoy - Qoxx + ^2^0 = 0, (20b) 

irot + iroy - tqxx + ^2^0 0, (20c) 

-isQt - isoy - soxx + L2S0 = 0, (20d) 

Lit = 0. (20e) 

Making use of (17), we rewrite (20e) in the following trilinear form 

(t>t4>xxt - 4>x4>tx4>tt + (Plt4>t + 4>x4>ttx4>t = 0. (21) 

The above trilinear equation ensures that the arbitrary manifold (f){x, y, t) should be parti- 
tioned as 

<t> = M^,y) + My,t), (22) 

where (f)i{x, y) and 02(?/, i) are arbitrary functions in the indicated variables. Making use of 
(22) in eqs. (18a) and (18b), one can show that eqs. (20a-20d) are consistent provided the 
submanifold ^2(2/, ^) can be split as 

02(y,t) = 02i(2/) + 022(t-2/), (23) 

Again, collecting the coefficients of (0°, 0°, 0°, (jp, (fP), we have only one equation 

L2t = 0. (24) 

Making use of (20a) for L2, (24) reduces to the form 

{Fut + Fuy)Fi + (Fu + Fiy)Fu = 0. (25) 

Equation (25) can be solved to obtain the form for Fi{y,t) as 

Fi{y,t)^Fi{t-y). (26) 
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Thus, the solutions of 2CLSRI can be written as 
S^'\x,y,t) 



where 



(01 (x, y) + 021 (y) + 022 (^ -y)y 



r 1 , 
J 2 <6i^ 



c(2). ,^ ^ V(022,-Fl(t-|/)^)0i.e^^ ^1 

^"^'^''^ (0l(x,y) + 02l(2/) + 022(t-Z/)) 



(0i(x,y) + 02i(|/) + 022(t-y))' 

201a;a; 

(01 (X, y) + 021 (l/) + 022 (t - y)) 



1 f {4>lyy + 4>21,yy) ~ ^01j;j;y 

2 V 01x 
(01y + 021,y) - ^01:r x 

1 ^01x2/ + 1 (01y + 021,y)^ + 



(27a) 



(27b) 



20L 

+ 022(t-y))' 

+ L2, (28) 



^2 01. 4 0L • ^''^ 

Thus, by choosing the arbitrary functions Fi{t—y), (f)i{x, y), 02i(?/) and (f)22{t—y) suitably, 
one can generate various kinds of localized solutions for the short waves S^^^ and S^'^^ while 
the longwave L does not support completely locahzed solutions. Prom (27a) and (27b), 
it is also obvious that the two physical fields S^^'^ and 5''-^-' have the same form except 
that their amplitudes are different and are driven by arbitrary functions Fi{t — y) and 
•\/022,t — -^i(^ — yY, respectively. It is also obvious that the 2CLSRI equation possesses an 
extra arbitrary function of space and time in comparison with its scalar counterpart [13]. 



IV. DROMION SOLUTIONS AND THEIR INTERACTIONS 

Now we choose specific forms of the arbitrary functions in (28) and (29) and obtain 
explicit exponentially localized dromion solutions and study their interactions. To generate 
a (1,1) dromion for the modes S'*^^-* and S^'^\ we choose the lower dimensional arbitrary 
functions of space and time, for example, as 
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Flit -y) ^ aisech{di{t - y) + ei) + gi 
(j)i{x, y) = a2tanh{b2X + C2y + 62) + g2 (30) 
hiiy) = ^3, '^22 - y) = a4anh{di{t - y) + 64) + 

Then, the corresponding exponentially localized solutions for and can be 

written as 



{aisech{di{t -y) + Ci) + ^1)^0262 sec /i(b2a; + Cs?/ + 63)^ 
(02 tanh(62a; + C2y + 62) + g2 + gs + 0,4 tanh((i4(^ - + 64) + 5(4)^ ' 

[a4(i4 sec h{di{t — y) + 64)^ — (oi sec h{di{t — y) + ei) + gii)^] x 

1^(2) ,2 _ "2&2SecM&2X + C2?/ + e2)^ 

(a2 tanh(62X + C2y + 62) + 6^2 + 5'3 + 04 tanh((i4(i - y) + 64) + ^4 
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The variable L takes the form 

2a|sec/i(62 + C22/ + 62)^ 



(a2 tanh(62a; + C2y + 62) +92 + 93 + 0,4 tanh((i4(t - y) + 64) + ^'4)2 
4a26| sec h{b2 + C2V + e^Y tanh(62 + + 62) 
(a2 tanh(62a: + C2y + 62) +92 + 93 + a4 tanh(d4(t - y) + 64) + 94) 

+ (32) 

A plot of the one dromion solution for the modes 5*^^^ and iS^^^ for the following parametric 
choice, «! = 1; 02 = 04 = 0.6; &2 = 1; C2 = 1; rfi = ^4 = 4; ei = 62 = 64 = 0; gi = 92 = 94 = 
0] 93 = 3(04(^4 > of) is shown in figs. 1(a) and 1(b). From the figures, it is clear that the 
dromions for the modes S^^^ and S'^^^ moving in the y-direction have different amplitudes 
and the amplitude of the dromions and hence the energy in a given mode depends on 
the parameter ai. We call such exponentially localized solutions driving S^^^ and 5*^^^ as 
"multimode dromions". Further, the above choice of lower dimensional arbitrary functions 
of space and time given by eq. (30) yields a line soliton for the long wave L as shown in fig. 
1(c). 

To generate a (2, l)-dromion for S'^^^ and S^'^\ we choose 

Fi = aisech{di{t — y) + ei) + 91 
01 = a2tanh{b2X + C2y + 62) + 

astanhibsx + Csy + 63) + 92 (33) 
021 = 93, 4>22 = a4tanh{d4{t - y) + 64) + 94 

so that the explicit solution can be written as 

^, N, .^[a2b2sech{b2X + C2y + e2y 
(ai sec h{di{t - y) + ei) + 9iy 

+0363 sec h(b3X + C3y + e^f] 

[02 tanh(62-« + C2y + 62) + 03 tanh(63X + csy + 63) + 92 + 93+ 
a4 tanh((i4(t — y) + 64) + 94]"^ 

[a4d4 sec h{d4{t — y) + 64)^ — (ai sec h{di{t — y) + ei) + yiY] x 

^ [«2^>2 sec h{b2X + C2y + 62)^ + 0363 sec h{b3X + Cay + 63)^] 
\S^^>\^ ^ . (34) 

[a2 tanh(62a; + C2y + 62) + 03 tanh(63X + c^y + 63) +92 + 93+ 



04 tanh(d4(i — y) + 64) + 94] 
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FIG. 2: Intensity profiles of the two dromion solution for the first mode (4a-4d) and second mode 
(4e-4h) at t = -6.0, -4.0, -1.0, 5.0. 

The plot of the (2,1) dromion solution for the modes S**-^-* and S**-^-* for the following 
parametric choice ai = 1; a2 = l;a3 = l;a4 = l;b2 = 63 = 1; C2 = 1;C3 = —l;di = — 
4; ei = 0; 62 = 63 = 0; 64 — gi — 92 — — 0; gz — 10(04(^4 > a\) at t=-6, -4, -1, 5 is shown 
in figs. (2a-2h). 

Prom the interaction of the dromions for the modes 5*^^^ and 5*^^) shown in figs. (2a-2h), 
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one observes that the two exponentially localized solutions with initial intensities Di and 
D2 {Di > D2) move along the diagonals in the (x-y) plane and exchange their intensities 
(energies) among themselves after interaction {Di < D2) thereby undergoing intramodal 
inelastic collision. It is also interesting to note that there is no exchange of energy between 
the two constituent modes and the energy contained in a given mode remains a constant. 

It should be mentioned that the choice of the lower dimensional arbitrary functions of 
space and time Fi{t — y), (f)i{x, y), (/>2(y) and 022(^ — y) determine the nature of the solutions 
admitted by 2CLSRI equation and their collision dynamics. For the choice of arbitrary 
functions given by eq.(30), one observes that the short waves are driven by exponentially 
localized solutions (dromions) and the energy contained in the first mode S^-^^ depends on 
Fi{t — yY while for the second mode S^^^ ^ it is governed by (022,t — — uY)- For the 
choice given by eq.(30) (with the parameters as in fig.4), the amplitude (energy) of the first 
mode 5"^^) is governed by the one dimensional soliton sech?A{t — y) while for the second 
mode S^'^\ it depends on the soliton 3sec/i^4(t — y). Thus, the choice given by eq. (30) 
launches two different energies in the modes 5*^^^ and 5*^^^ governed by the one dimensional 
solitons sec/i^4(t — y) and 3sec/i^4(t — y), respectively, and since the amplitude (energy) of 
the solitons docs not change during evolution, the energy contained in a mode remains a 
constant. Quantitatively, this is governed by the condition 

|5'(l)|2 ^ g2 

WWl^^i = constant (35) 

The above condition explains the existence of a firewall across the modes. This prohibition 
of energy across the modes by virtue of the existence of a firewall is valid only for the choice 
given by eq.(30), particularly if the short waves are to be driven by dromions. This behaviour 
in a vector (2+1) dimensional nonlinear pde is in sharp contrast to the Manakov model, a 
vector (1+1) nonlinear Schrodinger equation wherein the energy associated with the one 
dimensional solitons keeps flowing from one mode to the other. It should be mentioned that 
we report for the flrst time the identiflcation of exponentially localized solutions in a vector 
(2+1) dimensional nonlinear pde and their collision dynamics. 

From eqns. (27a) and (27b), one also observes that the sum of the squares of the short 
waves iS^^-* and 5"^^-* obeys the following equation 

' ' ' ' {(j)^{x,y) + (t>2i{y) + Mt-v)Y 
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(a) Dl D2 D2 Dl (b) 




FIG. 3: Time evolution of the composite mode S^^'"^^ at t= —6.0, 5.0 

where we call S^^''^^ as the "composite mode". Thus, we find that the composite mode 
is again driven by a two dromion solution (shown in figs.3((a),(b)) at t=— 6, 5). The 
intensity of the solution for the composite mode S^^'^^ is the sum of the constituent modes 
S^^^ and 5^^^ at every instant of time and one also observes a similar inelastic collision in 
the composite mode 5"*^^'^^ 

V. CONCLUSION 

In this paper, we have investigated the two component LSRI equation and shown that it 
admits Painleve property. We have then suitably exploited the truncated Painleve approach 
and generated multimode dromions for the short waves S^^^ and S^"^^ . The coUision dynamics 
of multimode dromions generated in the paper indicates that they suffer intramodal inelastic 
collision while the existence of a firewall prevents the fiow of energy from one mode to the 
other. It would be interesting to investigate the n-component LSRI equation from the 
perspective of localized solutions and their interaction. 
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